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' We examine the validity of the notion of the coherent state for pions and 

m ' 

0^ ' the quantum scattering effect in the final state of pion emission. When the 



O j' number of particles is large, the effect caused by the small but finite mass dif- 

43 



ference between the neutral and charged pions can add up substantially in the 



quantum evolution of an initially coherent state. As a result, the states with 
^ ■ quite different numbers of neutral or charged pions are essentially incoherent. 

The importance of the quantum scattering in the final-state isospin charge 
distribution of a disoriented chiral condensate (DCC) is investigated. We find 
that the scattering effect significantly reduces the spectacular Centauro and 
anti-Centauro events. The deformation of a charge distribution dP/df pre- 
dicted by the classical field theory is significant only for a DCC with the size 
of 10 fm or more. 
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I. INTRODUCTION 



Recently there has been a growth of interest in the subject of coherent pion radiation, 
in conjunction with the conjectured disoriented chiral condensate (DCC) It is sug- 

gested that the very high energy hadron-hadron or nucleus-nucleus collisions may create an 
extended region of space-time containing strong-interaction vacuum with a non-standard 
chiral orientation. In the language of the linear a-model the order parameter which usually 
points in the a direction is presumed to be in some other direction. In the "Baked Alaska" 
scenario 0, the disoriented vacuum is found in the expanding shell of hot collision debris. 
Indeed, the numerical simulations of the classical equation of motion show some evidence of 
growth of long wavelength modes when the expanding system is out of equilibrium, leading 
to the formation of large isospin correlated domains ||. When the hot shell hadronizes and 
breaks up, the disorientation is radiated away in Goldstone modes (pions) of a common 
isospin. The most natural way to connect the classical wave with a quantum state is the 
coherent state formalism in the limit of a large number of quanta. Such coherent pulses of 
semiclassical pion fields would lead to anomalously large fluctuations event-to-event in the 
ratio of neutral to charged pions produced. Assuming the all cartesian isospin directions are 
equally probable, one derives the probability (P) distribution 

- = — a) 

with the definition of neutral fraction / = N^o/N. 

In this paper, we examine the validity of the notion of the coherent state for pions, in 
particular, the quantum disturbance to coherence in the final state of pion emission. We 
find that the effect caused by the small but finite mass difference between the neutral and 
charged pions adds up substantially in the quantum evolution of an initially coherent state 
of pions. As a result, the states with quite different numbers of neutral or charged pions 
are essentially incoherent due to the different phase shifts! This phenomenon occurs only in 
the system where the number of isospin carrying particles (e.g., the pions) are large and the 
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evolution duration of the system is long compared to 0{l/m n ). Although the notion of the 
isospin coherence for a many-pion system can be entirely inadequate, the prediction of the 
neutral pion distribution P(f) of DCC in ([]]) is not altered by the phase shift caused only by 
the mass terms. We shall compute the modification of ([[]) due to the quantum rescattering 
effects when the interactions are included. 

II. ISOSPIN COHERENCE OF DCC 

The simplest quantum description of a DCC state is given in terms of the coherent 
state of the classical pion and a fields in the linear a model. When the DCC of a large 
isospin correlated domain is created, one considers an ideal limit where there is no spatial 
dependence of the DCC chiral orientation. Since we expect that only small isospin states are 
created, we consider for simplicity the 1 = state which is constructed by the superposition 
of the oriented DCC states denoted by \0,{p) over the isospin direction (9,(p) ||: 

|/ = 0) DC c = / dn(8,<p)\6,<p). (2) 

The DCC state relaxes to the true vacuum state by radiating the Goldstone modes. Note 
that the light particles emitted during this period are not the Goldstone modes associated 
with the true vacuum, but the disoriented vacua, which are in general linear combinations 
of 7i and a (they are "disoriented" Glodstone bosons). The classical a field evolves as well 
as the pion field, eventually the \I = 0)dcc state turns into the true vacuum (a constant o 
field) plus many outgoing low-energy pions. 

If one assumes that all pions radiated from a DCC during the relaxation end up in the 
same orbital state in the true vacuum, the state to which \I = 0)dcc relaxes can be written 
in a unique, simple form ]5],|6|]. The projection of this state onto an eigenstate of the pion 
number operator N = N + + AL + No is 

\I = 0; N) = , 1 {2ajal - ajaj)"/ 2 ^), (3) 
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where |0) is the true vacuum in which (a) = f n , and a± and a are the pion creation 
operators for and n°. The same neutral pion distribution P(f) as in ([!]) is obtained by 
expanding (^[) in the basis of definite charge states |Af )'s for / = N /N in the large N, N 
limit. Note that N + = AL = (N — Nq)/2 for \N$). As long as one accepts the assumption of 
a single orbital wavefunction for pions, the state \I — 0; N) is not subject to any final-state 
interaction corrections (up to an overall phase) since it is the only 1 = state made of A" 
pions. 

However, such a picture is oversimplified. It is not clear how all final pions in the 
true vacuum can retain the shape of the original wavefunction after a continuous emission 
of the time-varying "disoriented" Goldstone modes. These Goldstone modes emitted in a 
disoriented vacuum background field are different from those of the true vacuum in that 
they contain the ex's as well, which must turn into pions when the disoriented background 
changes into the true vacuum surrounding the DCC domain. As a result, more than one 
orbital wavefunction must be allowed for pions and many other isosinglet A" pion states can 
be constructed. Furthermore, the quantum scattering causes the transitions among those 
1 = states, and the coherence among different \N ) states may be lost. 

Even if all final pions end up in the same orbital configuration, the isospin breaking 
effects will break up the coherence among |Ao)'s. The isospin breaking caused by the 
electromagnetic interactions and the u- and d- quark mass difference is usually ignored in 
the DCC state and in its time development, which is a good approximation in most cases. 
However, it can be a problem when one is concerned with the coherence among states |A" ) 
with different A^'s in \I = 0;N). The ir^ and tt° actually oscillate with slightly different 
frequencies because of their electromagnetic mass difference Am, = m(7r ± ) — m(7r°) = 4.6 
MeV. Consider the single momentum mode k = 0. a+al|0) and the terms aJa^O) acquire a 
relative phase 8(t) = 2Am 7r t due to Am x , so that (up to an overall phase) 

\t; N) = . 1 {2a\al - e w <*> a\ a\) N ' 2 \0). (4) 

At t = 4 fm/c, 5(t) ~ 0.2. Although the phase for each 7r° pair is small, it accumulates with 
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Nq. The relative phase between the two charge states \Nq) and \Nq + ANq), differing in the 
number of neutral pions by AiVo, gets amplified by ANq times and becomes 5(t)AN . When 
N is let to infinity with a fixed Af = ANq/N, this phase goes to infinity and wipes out any 
coherence between the two states. For finite N and ANq, the initial coherence between \Nq) 
and \N + AiVo) is lost in a characteristic time scale r ~ 1/ (AN Am n ). For Af = 0.2 and 
N = 50, for instance, the isospin breaking effect rapidly accumulates in r ~ 4 fm/c, leaving 
only the I3 as a good quantum number. When the mutual coherence among different \Nq) 
states in (f|) is lost, and 7r° thereafter behave without the overall isosinglet constraint. 
However, it is interesting to note that the phase 8(t) does not affect the classical prediction 
(H) on the distribution P(f) |7j]. Our picture is very close to a purely classical description 
where the states with definite (classical) isospin orientation are produced incoherently, and 
the charge ratio N + /N, N_/N and Nq/N are identified with the isospin direction of the 
classical pion field. 



III. QUANTUM RESCATTERING EFFECTS 

We have argued that after the relaxation of the DCC into the true vacuum, the different 
charge states |iVo)'s should be treated as mutually incoherent. As we shall study how 
the final state quantum scattering changes the prediction in (|l|), we must start with a 
pion-number eigenstate [ AT ) instead of the isosinglet state in (||). The quantum scattering 
processes 7r°7r° <-> 7r + 7r~ cause the transition between \Nq) and \Nq ± 2), which can distort 
the distribution P(f). One needs to worry if there will be enough quantum scattering of 
pions to dilute the effect of interest in ([]]). Put more explicitly, suppose a DCC charge state 
initially consists of only iV 7r°'s, we would like to calculate how many 7r + 's or tt~'s end up 
in the detector after the quantum scatterings are included. 

The relevant quantity for determining the importance of the quantum scattering is the 
mean free path (Iq) of a physical pion in a cloud of emitted pions at the time (say, t = 0) 
when the DCC decays. If the size of the system (Rd) at the time is much smaller than l , 
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the pions essentially decouple from the system and propagate freely with little interactions. 
The opposite limit is Ru 3> Iq, in which case a pion will encounter many rescatterings before 
leaving the system and the measured value / (~ 1/3) has little trace of the original charge 
fraction Jq. In the general case, we have to conduct a quantitative computation. Since the 
pions from a DCC are soft in general, the mean free path of collision is much longer than 
a typical size of DCC. We shall make the approximation that the emitted pions stream out 
almost freely, i.e. that the collision does not slow down the diffusion of the pion cloud. In this 
approximation, there is a simple relation between the speed of diffusion (vd) and the average 
relative velocity (i> re i) for pion collisions. With this simplification (an underestimate), we 
shall show that the quantum corrections are characterized by a rather simple and intuitive 
factor exp[— 4R D /(3l )]. If the size of the diffusing cloud is a significant fraction of the 
mean free path, the quantum rescattering is non-negligible and will dilute the Centauro and 
anti-Centauro events 0. 

We are concerned with the evolution of the charge composition of the pion cloud with 
three charge states denoted by 0,+,—. For simplicity we assume that the pion density is 
spatially uniform inside the DCC and the typical momentum is determined by the uncer- 
tainty principle (p) ~ 1/Rd- For the DCC of an uniformly aligned orientation, this is a 
good estimate, but it is an underestimate for the DCC having domain structure. When 
the typical energy of pions is low, the relevant reaction is the two-body mr scattering. The 
neutral pion state, for instance, is depleted by the annihilation process 7r°7r° — > tt + tt~ and 
replenished by the reverse process 7r + 7r~ — > 7r°7r°. Under these assumptions, the time evo- 
lution of the numbers of three charge particles can be readily derived from the Boltzmann 
transport equation: 



dt 
dN. 



dt 
dN. 



dN, 



dt 



+ 








iV + iV_, 



(6) 



(7) 



(5) 
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where the cross section is defined to be that of 7r + 7r~ — > 7r°7r° (note that <j(tt + 7t~ — > 
7r°7r ) = icr(7r°7r° — > 7r + 7r~)). V^(£) is the volume of pion cloud enclosed by the outgoing 
pion front. (v TC \a n7T ) is the average value of the relative velocity times the cross section. 
Adding the above three equations gives the total number conservation dN/dt = where 
N = No(t) + N + (t) + N-(t) , which is the consequence of the fact that only two-body reactions 
are considered in collision terms in the Boltzmann equation. Since initially iV + (0) = AL(O), 
the evolution equations guarantee that this remains true at later times. In this case, we 
obtain a simple analytic solution. Define fo = N (0)/N and fit) = No(t)/N. The solution 
is 

(/o + l) + (3/o-l)e-^> 
m 3(/ + 1) - (3/o - l)e-Ht) ' W 

where Tit) determines the time dependence of charge composition due to the final-state 

scatterings, 

Ht) = 2N f Q (9) 

Before we estimate the Fit), some remarks on the general features of (H) are in order. 
Fit) is a positive number and e'^^ < 1. Generally / depends on J-(t) and the initial 
value f . It is easy to verify that when f < 1/3, / > f and when f > 1/3, / < f 
at later times (t > 0). In other words, if the system initially has fewer 7r°'s compared to 
7r + 's or 7r~'s, the quantum scattering tends to increase the number of 7r°'s; if more 7T°'s, 
the quantum scattering tends to decrease it. / = 1/3 is a special value, in which case the 
time dependence of / completely disappears and / keeps a constant / = / = 1/3. In the 
classical theory, f ranges from to 1 and the end points correspond to the Centauro and 
anti-Centauro events. These extreme situations do not occur in the quantum theory. In 
fact, / can never be equal to or 1 as long as J-{t) ^ 0. Since the state of all 7r° or all 7r =l= is 
degraded by the charge conversion processes, the quantum rescattering squeezes the allowed 
range for / from the both ends. The upper and lower limits on / are: 

1 + e~ T 1 — e~ T 

/max 7^ ' /min 7^ '. ~yc '■ (10) 

3 — e * 3 + e^ 
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The ranges for / < f min and / > / max are excluded by the quantum rescattering. If the 
rescattering is very strong, the factor e~ T would vanish, and / max and / m j n converge to a 
single point / = 1/3. Even for relatively weak scatterings of interest, the most spectacular 
Centauro and anti-Centauro events are diluted to some extent. 

We would like to estimate the integral in (^j for a realistic DCC. We assume that the 
volume V(t) of the pion cloud is spherical and diffuses with a constant velocity v& as 

V{t) = —{R D + v d tf. (0<t<oo) (11) 

The time dependence of the tttt cross section is negligibly weak as long as the typical momen- 
tum is in range of validity of the low energy theorem |J. We approximate it by the threshold 
value (by doing so we underestimate the effect). According to the soft-pion theorem, the 
model-independent low energy cross section is given by 

a ( n + n - ^ ir^ ) = ( 1 - ^ V , (12) 

whose threshold value is a n7T \ threshold = 9m£/1287r/^. The integration over t in flj^) can be 
done trivially, 

n(0)(v re ia n7T }R D 
JF(oo) = , (13) 

where n(0) = N/V(0) is the initial total pion density at the time when the DCC decays. 
Clearly, jF(oo) sensitively depends on the average relative velocity defined in a frame- 
independent way f rc i = y (pi ■ P2) 2 — rr^j E\E%. If the collision occurs frequently during 
diffusion, the diffusion velocity would be much smaller than the average pion velocity. In 
the other limit, if the collision rarely occurs, the diffusion takes place as pions stream freely. 
In the latter case is equal to (\p\)/E, and the average relative velocity is calculated by 
integrating i> re i over all relative directions: (v Te \) — ^(\p\)/E in the non-relativistic limit if 
all pions have approximately common Thus, we obtain 

Hoo) * (14) 
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where the mean free path at t = is Iq = 1/ (n(0)cr T7r ). The initial pion density is determined 
by the energy density of the DCC divided by the average pion energy, n(0) = m^f^/m n in 
the soft pion limit. Therefore the mean free path 1$ for the conversion 7r + 7r~ — ► 7r°7r° at the 
threshold energy is given by Iq ~ 1287r/^/9m^ ~ 27.7 fm, where /„- = 93 MeV. It is longer 
than the typical size of the DCC. As the cloud of quanta diffuses, the ratio of the mean free 
path to the cloud size increases even further. Therefore, for most DCC's, the diffusion of the 
pion cloud is closer to a free streaming than to a strongly interacting case. We thus justify 
the approximation in calculating the diffusion velocity va and the average pion velocity. For 
a DCC with a radius Rd = 5fm, jF(oo) is about 0.24. If the DCC consists of only 7r°'s, i.e. 
/o = 1, then /(oo) ~ 80%, i.e. 20% of original 7r°'s are converted into n^'sl 

We can also calculate the distribution dP/df of the final neutral charge fraction / = 
Nq{oo)/N. The classical field theory predicts dP/df = l/2y7o where f = N (0)/N is the 
initial value at the time when the DCC decays. The distribution in terms of / is given by 

dP _ dP df _ 1 16e _jr 

W~dfo'W~ 2^fW) [(l + 3 e -^)-3/(l- e -^)] 2 ' (15) 

where df /df is computed from (j|). The distribution dP/df tends to concentrate more in 
the range around / = 1/3 and fall off faster from one end to the other as / increases, 
making the distinction from non-DCC events a little less conspicuous. In Fig. 1 we have 
plotted dP/df for two values: Rd = 4 fm and 10 fm. The deformation is not too severe 
to significantly affect the experimental test of (|]) for Rd = 4 fm, but severe enough to 
make us worry for Rd = 10 fm. For values larger than Rd = 10 fm, the mean free path 
of collision becomes comparable or even shorter than the DCC size. Our assumption that 
(v re \) ~ 4v,i/3 breaks down and becomes slower than the average pion velocity. In other 
words, pions emitted stick together through final-state interactions. When this happens, the 
dilution factor !F(oo) rises quickly with Rd, making the experimental test of the classical 
theory prediction difficult. 
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Figure Captions 

Fig. 1 Comparison of the probability distributions of the neutral pion fraction / = 
iV"^o/iV tota i for different sizes of DCC 
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